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The correlated motion of electrons in multi-orbital metallic ferromagnets is in- 
vestigated in terms of a realistic Hubbard model with AA-fold orbital degener- 
acy and arbitrary intra- and inter-orbital Coulomb interactions U and J using a 
Goldstone-mode-preserving non-perturbative scheme. An effective quantum parame- 
ter 'h' = ]^jj^^^Z\^jyi. is obtained which determines, in analogy with 1/5 for quantum 
spin systems and 1 /N for the A^-orbital Hubbard model, the strength of correlation- 
induced quantum corrections to magnetic excitations. The rapid suppression of 
this quantum parameter with Hund's coupling J, especially for large AA, provides 
fundamental insight into the phenomenon of strong stabilization of metallic ferro- 
magnetism by orbital degeneracy and Hund's coupling. This approach is illustrated 
for the case of ferromagnetic iron and the half metallic Heusler alloy Co2MnSi. For 
realistic values for iron, the calculated spin stiffness and Curie temperature values 
obtained are in quantitative agreement with measurements. Significantly, the con- 
tribution of long wavelength modes is shown to yield a nearly ~ 25% reduction 
in the calculated Curie temperature. Finally, an outline is presented for extending 
the approach to generic multi-band metallic ferromagnets including realistic band- 
structure features of non-degenerate orbitals and inter-orbital hopping as obtained 
from LDA calculations. 



PACS numbers: 71.10.Fd,75.10.Lp,75.30.Ds,75.40.Gb 



2 



I. INTRODUCTION 

Dramatic improvements in experimental techniques such as Angle Resolved Photoe- 
mission Spectroscopy (ARPES)^ have led to important insight into the origin and role 
of correlation effects in itinerant ferromagnets such as iron, highlighting the coupling of 
electrons with magnons as playing a major role in the electron self energy renormalization 
and scattering rates. Correlation effects also play an important role in the observed zone 
boundary magnon softening and damping in ultrathin films of iron, as observed in recent 
Spin Resolved Electron Energy Loss Spectroscopy (SPEELS) experiments,- where the 
zone boundary magnon energies are much lower than those predicted within the random 
phase approximation (RPA), and the magnon energies were observed to depend non- 
monotonically on the film thickness.- Further evidence of correlation effects is provided 
by ab-initio band structure calculations of half-metallic Heusler alloys showing emergence 
of non-quasiparticle (NQP) minority-spin states near the Fermi energy at finite tempera- 
tures, which has been suggested to be responsible for the strong suppression of tunneling 
magnetoresistance (TMR) ratio with temperature in Co2MnSi-based magnetic tunneling 
junctions (MTJ) observed in tunneling conductance measurements.- These experiments 
conclusively highlight the importance of incorporating electron-magnon coupling effects 
in the correlated electron spin dynamics in metallic ferromagnets. 

A considerable amount of work has been devoted to understanding the electronic band 
structure of these metallic ferromagnets,- starting from the Local Spin Density Approxi- 
mation (LSDA) within the Density Functional Theory (DFT),- which accounted for corre- 
lations only in a limited way. The development of several extensions such as the LDA-I-U, 
LDA++, and LDA+DMFT,^-^ has led to considerable progress in incorporating correla- 
tion effects in realistic band structure calculations. However, here the correlation term 
is incorporated either at the mean-field level or within a local self-energy approximation 
which neglects the momentum dependence. These methods therefore cannot be used to 
directly address spin wave excitations as they do not explicitly preserve the spin rotation 
symmetry, for which vertex corrections must also be included systematically,— and also 
predict much higher Curie temperatures^ than observed experimentally due to neglect of 
long wavelength spin fluctuation modes. 

Since metallic ferromagnets are characterized by intermediate to strong correlations, a 
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proper description of spin waves must incorporate correlation effects non-perturbatively 
and simultaneously preserve the Goldstone mode. There are mainly two theoretical ap- 
proaches for studying spin- wave excitations in itinerant ferromagnets — the random phase 
approximation (RPA)/^"^^ and mapping to an equivalent Heisenberg model of localized 
spins by using the magnetic force theorem and its generalization to compute the exchange 
interaction parameters.— Due to neglect of correlation effects, the RPA is well known 
to overestimate the spin stiffness, magnon energies, and stability of the ferromagnetic 
state.— On the other hand, mapping to an effective Heisenberg model does not capture 
typically itinerant features such as zero-temperature magnon damping. The adiabatic 
approximation has been used to investigate spin dynamics of ultra-thin films, but it 
has been pointed out that this approach breaks down for large wave-vector modes.— 
The localized spin model has also proved unsatisfactory in explaining the doping de- 
pendence of the anomalous softening and damping of zone boundary spin wave modes in 
the OMR manganites.— Although Linear Response Density Functional Theory (LRDFT)- 
based studies of spin dynamics of irort^'^ and the Heusler alloys^° account for damping of 
high-energy magnon modes due to decay into Stoner excitations, it has been pointed out 
that spin-charge coupling in a band ferromagnet results in significant magnon damping 
for modes lying even within the Stoner gap.— 

In this situation, it will be useful to have a scheme which could incorporate features 
of the realistic electronic band structure and simultaneously take into account the most 
important correlation effects within a non-perturbative and Goldstone mode preserving 
approach. Such a scheme would also be useful from the technological point of view, 
since many half-metallic ferromagnets such as the Heusler alloys are being intensively 
investigated due to their potential applications in the spintronics industry. 

Recently, correlation effects in metallic ferromagnets have been investigated using a 
non-perturbative, inverse- degeneracy based expansion scheme in which self energy and 
vertex corrections are included systematically so that the spin rotation symmetry and 
the Goldstone mode are explicitly preserved order by order.—"— Within this approach, 
the enhancement of ferromagnetism due to suppression of correlation-induced quantum 
corrections was found to be strongly dependent on several electronic band features such 
as lattice-type, dimensionality, t'-induced DOS asymmetry, and band filling. The magnon 
self-energy was also investigated in the context of zone boundary magnon softening and 
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magnon damping observed in manganites and ultra-thin transition-metal films.— i^si^^ 

The A/'-orbital Hubbard model considered in the above works involved the orbitally- 
symmetric case with identical intra-orbital and inter-orbital Coulomb interactions. In 
this paper, we shall extend this spin-rotationally symmetric Goldstone-mode preserving 
approach to the general orbitally-asymmetric case with arbitrary intra- and inter-orbital 
Coulomb interactions. This provides an important extension of our recent work on quan- 
tum corrections in the doubly-degenerate {JV = 2) Hubbard model with intra-orbital 
Coulomb interaction U and inter-orbital Hund's exchange J.— We will derive an effective 
quantum parameter 'h' = {U^ + (A/" — 1) J^) /{U+ (A/" — 1) J)^ which determines, in analogy 
with l/A/" for the orbitally-symmetric Hubbard model and l/S* for quantum spin systems, 
the strength of the correlation-induced quantum corrections in a realistic multi-orbital 
metallic ferro magnet. 

As an illustration of this spin-rotationally invariant approach for investigating correla- 
tion effects in multi-band metallic ferromagnets, we will calculate electronic and magnetic 
properties in a five-orbital model with realistic parameters corresponding to ferromagnetic 
iron. The inclusion of long wavelength spin-fluctuation modes in the finite-temperature 
spin dynamics is distinct advantage of our approach, and it provides a quantitative mea- 
sure of the Curie-temperature overestimate in local self-energy based calculations (for 
example, the LDA+DMFT) which neglect contribution of long wavelength modes. We 
will also outline an extension to a realistic multi-band model including non-degenerate 
orbitals and inter-orbital hopping as typically obtained in a LDA calculation. 

The outline of this paper is as follows. After introducing the orbitally-asymmetric 
A/'-orbital Hubbard model in Section II, the transverse spin fluctuation propagator is 
studied in section HI, and the effective quantum parameter is derived from the first order 
quantum correction diagrams for the irreducible particle-hole propagator obtained here. 
Our approach is illustrated with two applications corresponding to ferromagnetic iron 
in Section IV and the Heusler alloy Co2MnSi in Section V. The extension to realistic 
multi-band model including non-degenerate orbitals and inter-orbital hopping is outlined 
in Section VI, and our conclusions are presented in section VII. 
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II. AA-ORBITAL HUBBARD MODEL WITH HUND'S COUPLING 

We consider the following orbitally asymmetric A/'-orbital Hubbard model with arbi- 
trary intra-orbital (U) and inter-orbital (J) Coulomb interactions: 



where /i and z/ refer to the TV degenerate orbitals at each lattice site i, and Sj^ = 
ip\^{(T /2)ipi^ is the local electron spin operator for the /i orbital in terms of the corre- 
sponding fermion operators ip\^ = {a]^^ a]^|^) and the Pauli matrices cr. For J = U, the 
model reduces to the orbitally symmetric case considered earlier.— The inter-orbital den- 
sity interaction term Vni^nii, is not included here as this charge term has no leading order 
effect on magnetism and only weak effects on spin dynamics when quantum corrections 
are included away from the onset of staggered orbital ordering. The role of orbital fluctua- 
tions on spin dynamics due to this inter-orbital density interaction term has been studied 
recently in the context of the observed zone-boundary anomalies in manganites.— 

The continuous spin rotation symmetry of the above Hamiltonian implies the existence 
of Goldstone modes in the spontaneously broken symmetry state. In the following we 
will present a non-perturbative scheme in which the correlation-induced self energy and 
vertex corrections are incorporated systematically so that the Goldstone mode is explicitly 
preserved order by order. 



We assume a ferromagnetic ground state with magnetization in the z direction and 
examine transverse spin fluctuations representing both collective (spin-wave) and single- 
particle (Stoner) excitations. We consider the time-ordered transverse spin-fluctuation 
propagator in this broken-symmetry state: 



where /i and v are any of the M orbital indices, and the fermion-spin lowering and spin- 
raising operators Sf = ip\ Aa^ /2)%lji^. 




(1) 



III. TRANSVERSE SPIN FLUCTUATIONS 



(2) 
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In terms of the irreducible particle-hole propagator 0^i^(q, ou), the spin-fluctuation prop- 
agator can be written exactly as: 

X/Ii/^(q>^) = 0/^1^ (q>^) + 0MM'(q>^)^MvXMV'(q>^)> (3) 

where the interaction term U/^i, — U ior /j, — i/, and U/^i, — J ior /i u, and summation 
over repeated indices is implied. It is physically relevant to consider the total transverse 
spin-fluctuation propagator: 

as it measures the response to external probes such as magnetic field or neutron magnetic 
moment which couple equally to the electron moment for all orbitals. It is particularly 
convenient to solve the coupled equations then to obtain: 

where the interaction term C/+ — U + {J\f — 1) J, and the total irreducible particle-hole 

propagator: 

0(q,w) = Xl0a,.(q,i^)- (6) 

In analogy with the 1/jV expansion for the orbitally symmetric j\A-orbital Hubbard 
model, we consider a systematic expansion: 

= 0(0) + 0(1) + 0(2) + . . . (7) 

for the irreducible propagator 0(q, u) in terms of fluctuations. The first term 0(°) is simply 
the bare particle- hole propagator, whereas the higher-order terms 0(^), 0(2^ etc. represent 
correlation-induced quantum corrections involving self energy and vertex corrections. 

A. Random Phase Approximation 



Retaining only the zeroth-order term 0( ) in the expansion yields the random phase 
approximation, amounting to a "classical-level" description of non-interacting spin- 
fluctuation modes. As the hopping term is diagonal in orbital indices, the zeroth-order 
term involves only the intra-orbital contribution: 

0W(q,(^) = xo(q,^) = Yl -J+ -' (S) 
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where the Hartree-Fock level band energies e^ = e^^ — a A involve the exchange splitting: 

2A = [f/+ (A/"- l)J]m (9) 

between the two spin bands. The superscripts + (— ) refer to particle (hole) states above 
(below) the Fermi energy ep- Here the magnetization m = 2(S'f^) is identical for all Af 
orbitals in the orbitally degenerate ferromagnetic state. For the saturated ferromagnet, 
the magnetization m is equal to the particle density n for each orbital. 

Due to orbital degeneracy, there are only two independent cases of interest correspond- 
ing to jj, = u and [i v. Following Eq. (|3]), the two corresponding coupled equations at 
the RPA level are: 

Xaa =X0 + X0f/Xaa +(A^-1)X0JX^+, (10) 

Xjt = XoJx~a + Xof/X^+ + (A/" - 2)xo^X^+, (11) 
solving which, we obtain 

where the two interaction terms above are = U + (A/" — l)J and U~ = U — J. 
The propagators involve linear combinations of in-phase and out-of-phase modes with 
respect to the orbitals, representing gapless (acoustic) and gapped (optical) branches, 
respectively. The in-phase mode with effective interaction corresponds to the usual 
Goldstone mode (acoustic branch), while the out-of-phase mode with effective interaction 
yields gapped excitations (optical branch).— 



B. Quantum corrections and effective quantum parameter 

Diagrammatic contributions to the first order quantum correction 0^^^ for the orbitally 
asymmetric Hamiltonian (1) are shown in Fig. [H Structurally, they are similar to the 
0(1/^/) diagrams in the orbitally symmetric case,^^ but the different orbital components 
with appropriate interaction terms are now considered separately. Diagrams (a) and (d) 
represent corrections to the irreducible propagator due to self-energy corrections, whereas 
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FIG. 1: First-order quantum corrections to the irreducible particle-hole propagator (j){ci,uj) for 
the realistic AA-orbital model with arbitrary intra-orbital and inter-orbital Coulomb interactions 
U and J. 

diagrams (b) and (c) represent vertex corrections. The corresponding expressions are 
obtained as:— 



E / ^[iu' + i^^-l)J')x-JiQ,^) + 2i^^-l)UJx-,:{Q,n) 

1 



*^k'-q+Q 



e^, + u — fl — ir] 



(14) 



2E / ^^{ur-JiQ,^) + i^f-^)JK^iQ,m 
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el„ + Q - ir] 



(15) 



= E / ^{[f/' + (A/'-i)^']r;.+(Q,^) + [2(A/'-i)Jt/ 

Q 

+ (Ar-l)(Ar-2)j2]r-;(Q,(])} 
1 
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+ Q - ir] 
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(16) 
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1 



4" + - ^ - 



(17) 



where the kernels and T for the spin-wave propagators are defined in terms of Xc 
and XaB^ from Eqs. fll2fl3p are obtained as: 



[r-+ 



1 _ [Xaa ]rPA - XO _ 1 

. _ [Xap]KPK _ 1 r U 



?7+ ^(Ar-l)f/- 



1 - f/+Xo 1 - U~Xo 



(19) 



l-f/+Xo l-f/-Xo. 

To demonstrate the exact cancellation and hence the Goldstone mode for g = 0, we 
note that the boson term (quantity in braces) in Eq. (fT5|) for 0*^^^ can be expressed as: 



UT-^ + {M- l)JT-; = + (A/- - 1) J2]x-+ 

+ [2{M-l)UJ+{M-l){M-2)j']x~J]/x^ 



(20) 



which is identical to the boson term in Eq. f|T^ for The above identity is shown in 
the Appendix. Similarly, using Eqs. ( JT8l) and ( IT9l) . the kernels in the boson term of Eq. 
( |T6l) for 0*^'^-' can be written in terms of x ^ ^"^^ Xo- Using the above substitutions in Eqs. 



10 



([T1-(II7D, and with ej;" 




2A for g = 0, we obtain: 




-2 



+ 



+ 



X 



+ (AT - l)j')x-: + (2(Ar - 1)UJ + (AT - l)(Ar - 2)j')x-,n 
{{U' + (AT - 1) J^)x-+ + (2(Ar - l)f/J + (AT - l)(Ar - 2)j')x-,n 
{{U' + (A/- - 1) J^)(x^+ - xo) + (2(A/' - 1) Jt/ + (AT - l)(Ar - 2) J2)xj+} 
{(f/' + (A/'-l)^')Xo}], 



(21) 



which yields identically vanishing contribution for each spin-fluctuation mode Q. We 
note that this mode-by-mode exact cancellation is quite independent of the spectral dis- 
tribution of the spin-fluctuation spectrum between collective spin-wave and particle-hole 
Stoner excitations. Furthermore, the cancellation holds for all u, indicating no spin- wave 
amplitude renormalization, as expected for the saturated ferromagnet in which there are 
no quantum corrections to magnetization. 

We shall now obtain an effective quantum parameter which approximately determines 
the strength of the quantum corrections obtained above. Due to the uncorrelated nature 
of the inter-orbital spin fluctuations {{S^S^^) = 0), the contribution of the inter-orbital 
propagator in Eq. ( 11^ is much smaller than the contribution of the intra-orbital 
propagator Xaa^ and hence only the contribution from the orbitally diagonal term [U"^ + 
(A/" — ^)J'^]Xaa essentially survives, leaving an overall factor U'^ + (A/" — 1)J^ on carrying 
out the Q, n integration. Comparing with the corresponding factor {U + (A/" — 1)^/)^ 
obtained for the equivalent single-orbital case (with identical exchange splitting) yields 
an overall relative factor of [f/^ + {Af — 1)J'^]/[U + {Af — l)J]'^- Similarly, the quantum 
corrections (p^^^ and (p^'^^ also yield the same overall factor. Thus the total first- 
order quantum correction approximately bears an overall relative factor of [f/^ + (A/" — 
1)J^]/[U + (A/" — 1) J]^ compared to the equivalent single-orbital case. This relative factor 
thus plays the role of an effective quantum parameter which determines the strength 
of the correlation-induced quantum corrections in a multi-band metallic ferromagnet. 
This quantum parameter is exact in the orbitally independent limit J/U ^ where it 
approaches 1, and in the orbitally symmetric limit J/U — i- 1 where it approaches l/A/". 
Also, the quantum parameter falls rapidly with Hund's coupling J, especially for large Af, 
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highlighting essentially the role of orbital degeneracy and Hund's coupling in stabilizing 
metallic ferromagnetism by suppressing the quantum corrections. 



C. Spin stiffness 

As the Goldstone mode is explicitly preserved in our approach, it allows investigation 
of correlation effects on spin stiffness and hence on the ferromagnetic stability with respect 
to long wavelength fluctuations. Here, we will evaluate the first-order quantum correction 
to spin stiffness exactly and then compare its J dependence with that of the effective 
quantum parameter obtained above. This quantitative comparison will clearly show the 
usefulness of the effective quantum parameter. 

First-order quantum corrections to spin stiffness are derived by expanding ^'■^''(q) for 
small q as in the two-orbital case.— There is no quantum correction to the delocalization 
contribution (V^ek) in the spin- stiffness constant; only the exchange contribution in the 
spin stiffness is renormalized by the surviving second-order terms in 5 = Ck — Ck-q, and 
we obtain for the first-order quantum correction to spin stiffness: 

D« =2A(f/+)0«/?' 
1 U+ 

d (2A)3 ^ J 2m 
Xo(Q,^^) 

^ f/eff(Q,^) 




+ [U^ + {M-i)J^] 

where the effective interactions: 



(22) 



u:^ = [u' + (AT - i)j']x-j + [2{Ar - i)uj + {Ar-i){Ar- 2) j^ix^^, 

u:^ = u:^-{u' + {M-i)j')xo. (23) 

The calculated quantum correction to spin stiffness, normalized so that it equals 1 for 
J/U = 0, is plotted in Fig. [2] as a function of J/U for A/" = 5 orbitals. Here we have 
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FIG. 2: The normalized first order quantum correction D^^\j) / D^^\0) to spin stiffness calcu- 
lated from Eq. (I22p compared to the effective quantum parameter [1 + {J\f — 1)( J/[/)^]/[l + 
(A/" - 1)( J/C/)]2 as a function of J/U for = 5 orbitals with fixed U + {M-1)J = 1.5W = ISt. 

considered the sc lattice, band filling n = 0.3, f = 0.25, and fixed U + {Af — 1)J = 1.5W. 
Also shown for comparison is the quantum parameter 'h' obtained above, which is seen to 
be exact in the orbitally independent {J/U — )■ 0) and the orbitally symmetric {J/U — )• 1) 
limits, and it remains close to the calculated corrections even in the intermediate region, 
as expected from the uncorrelated nature of inter-orbital spin fluctuations.— 

Including the quantum correction, the renormalized spin stiffness is then obtained as: 

D = - , (24) 

where the bare-level (RPA) spin stiffness: 

(25) 



= - 
d 



_2' ' 2A 

involves two characteristic terms representing delocalization energy loss and exchange 
energy gain upon spin twisting. 

IV. APPLICATION TO IRON 

In the previous section we obtained an effective quantum parameter iii 
terms of the physically important parameters A/", U, J for a multi-orbital band ferromag- 
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net. This quantum parameter determines, in analogy with l/A/" for the generahzed Hub- 
bard model and 1/5* for quantum spin systems, the strength of quantum corrections to 
magnetic excitation energies. The quantum parameter is strongly suppressed by Hund's 
coupling J, and rapidly approaches the limiting value of l/A/", especially for large M . 

This suggests that quantum corrections in a realistic multi-band ferromagnet with 
arbitrary A/", f/, J can be conveniently investigated in terms of the orbitally symmetric 
Hubbard model: 

H = -t «Iaa%-fTa + ^ X^(«lta«ita'^4/3«4/3 + «Ita"*t/3«4/^ (26) 

{ij),a,a i,a,P 

and the inverse-degeneracy expansion scheme,— with an effective number of orbitals 
such that = ^jj^lj^Zi)j)2 ■ III this section, we will apply this approach specifically to the 
case A/" = 5 with realistic parameters corresponding to ferromagnetic Fe. Accordingly, we 
consider a bcc lattice with dispersion: 

ek = %tcos{kxa)cos{kya)cos{kza) + 2t' cos(2fc^a), (27) 

where t and t' are the nearest and next-nearest neighbor hoppings, which yields bandwidth 
W = 16t (for t'/t < 0.5). We have taken the lattice parameter 2a = 2.87A, t'/t = 0.5, 
t = 0.2 eV, J/U ~ 1/4, and the interaction strength U of the order of bandwidth W as 
appropriate for a strongly correlated system. This yields U c:^ W = 16t = 3.2 eV and 
J = 0.8 eV, which are close to the parameter values considered in band structure^ and 
constrained LDA^ studies for iron. With these parameters, and N' = 5 corresponding 
to the five 3d orbitals in iron, the quantum parameter 'h' = ^u^lj^Zi)j)2 — |) which 
corresponds to ~ 3 within the A^-orbital Hubbard model. The choice of t'/t = 0.5 is 
motivated by the fact that in a bcc lattice, the second neighbors are only about 15% more 
distant than the nearest neighbors, and the hopping term tij depends on the inter-site 
distance as t^j oc 1/rf-,^^ which yields t'/t ^ 0.5. 



A. Spin stiffness 

Fig. [3] shows the renormalized spin stiffness for the bcc lattice for different number of 
orbitals A^. The spin stiffness is negative for = 1 and rapidly becomes positive with in- 
creasing A^, indicating the strong role of orbital degeneracy in stabilizing ferromagnetism. 
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FIG. 3: The renormalized spin stiffness as a function of band filling for different effective number 
of orbitals N, evaluated for the bcc lattice with realistic bandwidth W = 16t = 3.2 eV, Coulomb 
interaction energy U = W = 3.2 eV, and lattice parameter 2a = 2.87A for Fe. The measured 
value for Fe is 280 meVA^. 

The optimum filling occurs near n = 0.5 for finite N. With realistic parameters taken 
for Fe as given above, the calculated values of the renormalized stiffness are close to the 
measured value of 280 meVA^ for iron obtained from neutron scattering studies.— This 
indicates that a simple multi-band Hubbard model with an effective number of orbitals 
to incorporate the quantum corrections provides a quantitative description of the spin 
stiffness for strongly correlated metallic ferromagnets such as Fe. 



B. Renormalized magnon dispersion 

The renormalized magnon energy Wq for mode q was obtained from the pole condition 
[1 — f/Re0(q, — Wq) = 0] in the total spin fluctuation propagator: 

X q,W =l 77-7 7 28 

1 - f/0(q,w) 

where the irreducible particle-hole propagator: 

= 0W + (29) 

up to first order in 1/A^. The numerical evaluation of the quantum correction 0*^^^ by 
integrating over the intermediate (Q, fl) states has been discussed earlier.— Both collective 



15 



0.8 



0.6 



0.7 



U=3.2eV 
n=0.5 
t=0.2 eV 
t'=0.5t 



BCC 



{N=5) 

COq {N=3) » 

Mq{N=1) e 




0.5 



> 



3 







r 



X 



M 



R 



r 



FIG. 4: Renormalized magnon energies for the bcc lattice for different number of orbitals N, 
showing a rapid crossover from negative-energy to positive-energy long wavelength modes as 
N is increased from one, showing the strong role of orbital degeneracy in stabilizing metallic 
ferromagnetism. 

and Stoner excitations are included. While the bare particle-hole propagator (/)^°)(q, w) 
remains real in the relevant tu range, the quantum correction 0*^^^(q, co) is complex for any 
finite u due to the couphng with charge fluctuations, resulting in flnite zero temperature 
magnon damping.— 

Fig. m shows the renormalized magnon energy dispersion for the bcc lattice for different 
A^. The magnon energy exhibits a Goldstone mode at both F and M, as they are equivalent 
points in our extended Brillouin zone which extends from —Tr/a to n/a in each direction, 
whereas the bcc lattice parameter is 2a. We flnd that for = 3 and 5, the magnon energy 
renormalization is nearly momentum independent in the F-X, X-M, and M-R directions. 
However, near (7r/2, 7r/2, 7r/2) between F-R, the magnon energy is softened relatively more 
strongly. 



We now evaluate the renormalized electronic density of states (DOS) by incorporating 
the self energy correction due to electron-magnon coupling. A J,-spin particle (energy 
^ic > ^f) can decay into a magnon and a t-spin particle, resulting in considerable J,- 



C. Density of states and quasiparticle dispersion 
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spin spectral weight transfer to just above the Fermi energy. Since we are considering 
a saturated ferromagnet with no J,-spin density, an '|"-spin particle cannot decay into a 
4-spin particle due to spin conservation, and hence there are no quantum corrections at 
T = to the t-spin DOS. 

The J,-spin self-energy J^^ is calculated within an approximate resummation procedure 
which incorporates particle-particle correlations:— 

where eI^"* = 2A = mU is the HF-level self energy, and the first order self energy: 

= Imf/^y . (31) 

Q ^ "Q ^k-Q ^ 

Fig. [5] shows the renormalized density of states for =3 and 5 orbitals. The J,-spin 
band HF DOS is also shown for comparison. The 4,-spin DOS is seen to be renormal- 
ized considerably with significant band narrowing, band shift, as well as spectral- weight 
transfer from the bare band to just above the Fermi energy. 

Fig. |6] shows the energy momentum dispersion of ^-spin electrons in terms of an 
intensity plot of the renormalized spectral function A]ii{uj) obtained from the Green's 
function: 

Go (k,w) - S;(k,a;) 

The bare (HF) dispersion for the minority spin is also plotted for comparison. The renor- 
malization is seen to be especially strong near the X and R points, indicating significant 
band flattening and mass renormalization, and also strong non-quasi-particle character of 
the low-energy minority-spin states near the Fermi energy, as further discussed below. 

An important aspect in Fig. [5] is the emergence of new non-quasiparticle (NQP) 
states just above the Fermi energy at K, corresponding to strongly incoherent spectral 
function near the X and R points in Fig. [6l These NQP states are important in view of 
recent tunneling conductance measurements on Heusler alloy-based Magnetic Tunneling 
Junctions, where a strong suppression of spin polarization is observed with temperature.- 
This suppression has been attributed to the emergence of NQP states at and below the 
Fermi energy due to correlation effects arising from the electron-magnon coupling. 
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FIG. 5: The spin-resolved density of states for the bcc lattice at band filling n = 0.5, showing the 
transfer of ^-spin spectral weight due to correlation effects and the emergence of non-quasiparticle 
(NQP) states just above the Fermi energy. 

D. Finite temperature spin dynamics 

With correlation effects incorporated, the irreducible particle-hole propagator 0(q, u) 
yields renormalized magnon energies as obtained above, and also the effective spin cou- 
plings Jij = U'^(f)ij in an equivalent Heisenberg model after integrating out the fermionic 
degrees of freedom. At finite temperature, these renormalized magnon energies will de- 
termine magnetization reduction resulting from the electron-magnon coupling self energy 
which is explicitly of order 1/A^ as in Eq. ( |3T1) . This is also similar to the low-temperature 
relative magnetization reduction of order 1/S in the spin-S" Heisenberg ferromagnet. In 
vew of this equivalence with an effective spin-S* Heisenberg model with = 2S, we will 
therefore calculate the finite temperature magnetization from the self-consistent Callen 
equation for the spin-S" Heisenberg model:— 

^•^^/^ ~ (1 + $)25+l _ $25+1 ' ^^^^ 

where the magnon amplitude: 

$ = y -J^ , (34) 

Q 
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FIG. 6: (color online) Intensity plot of the renormalized 4,-spin spectral function (bottom) along 
different symmetry directions in the Brillouin zone for the bcc lattice (with = 5). Also shown 
for comparison is the 4,-spin bare quasiparticle dispersion (top). 



in terms of the thermally renormalized magnon energies: 

= ^ci{Sz)t/ {Sz)o- 



(35) 



Solving the above three coupled equations self-consistently with S = N/2 yields the 
temperature dependence of the magnetization. While the finite-temperature electronic 
density and magnetization corrections, which renormalize the u term in the magnon 
propagator, are incorporated in this calculation through the thermal renormalization of 
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FIG. 7: The RPA- level magnon DOS for the bcc lattice, showing a peak at the dominant mode 
energy fipp^ = 0.26 eV. 
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FIG. 8: Including long wavelength modes in the self-consistent Callen scheme for magnetization 
as a function of temperature results in a nearly 25% suppression of Tc, as seen by compari- 
son with the dominant mode approximation (DMA), highlighting the significant role of long 
wavelength spin fluctuations in determining the Curie temperature. 

magnon energies, it neglects intrinsic finite-temperature renormalization effects on the 
spin couplings. 

Instead of performing the computationally intensive Q summation in Eq. ( 134|) exactly, 
it is convenient to break it into two parts, corresponding to long- and short-wavelength 
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contributions: 

^=y^ — = y^ — — , (36) 

Q ^ Q<A ^ ^ Q>A ^ ^ 

where A represents a momentum space cut-off for the long wavelength modes. The small-Q 
contribution can be evaluated explicitly by integration, whereas the large-Q contribution 
can be approximated by a single term corresponding to the "dominant mode energy", 
since the magnon DOS exhibits a pronounced peak at this energy corresponding to short 
wavelength modes. Thus, $ = $lw + '^'dma, corresponding to the long wavelength con- 
tribution and the dominant mode approximation. The evaluation of $lw is discussed 
below. 

For small-Q modes, the magnon energy f^q = DQ^ in terms of the thermally renor- 
malized spin stiffness: 

b = D{S,)t/{S,)o (37) 

at temperature T, where D incorporates the zero-temperature quantum corrections dis- 
cussed above. Using the transformation x = DQ'^/k^T, we obtain: 

where Xc = DA^/k-gT corresponds to the momentum cut-off A. For low temperatures 
{ksT <C -DA^), the cutoff Xc ^ 1, and the x-integral is nearly temperature independent, 
which yields the well-known Bloch T^/^ law for the decay in magnetization. 

While including the long wavelength contribution, it should be noted that apart 
from (0,0,0), there are twelve other points in the extended Brillouin zone which cor- 
respond to the Goldstone mode. These points are located at (ivr, ivr, 0), (0, ±7r, ±7r), 
and (±7r,0, ±7r), and each of these points is associated with a quadrant of a sphere in 
momentum space. Therefore, there are effectively four points in the Brillouin zone around 
which Qq = DQ^, and hence the factor 4 in Eq. fl38|) . 

Fig. IH] shows the finite temperature magnetization evaluated by solving the three cou- 
pled equations ( 133115^ with 5 = 3/2 and the renormalized magnon energies obtained in 
subsection B. The results are shown for both the dominant mode approximation (DMA) 
and after including the long wavelength (LW) contribution. The momentum cut-off for 
long wavelength modes was fixed at Aa = 1. The Curie temperature obtained within the 
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FIG. 9: Lattice structure of Co2MnSi. 

DMA is approximately 1350 K. Including the long wavelength contribution reduces the 
Curie temperature to 1020 K. An important result of our analysis is that long wavelength 
modes yield a significant (nearly 25%) reduction in the Curie temperature. This compari- 
son provides a useful quantitative measure of the overestimation of the Curie temperature 
in calculations which neglect long wavelength modes such as in DMFT studies for iron.- 

V. APPLICATION TO CoaMnSi 

As another application, we consider the Heusler alloy Co2MnSi, which, owing to its 
half-metallicity and high Tc has attracted considerable attention in recent years due to 
possible applications in fabricating spintronics-based devices.— i^^"— The key feature being 
exploited is the 100% spin polarization arising due to the gap in the minority-spin DOS 
at the Fermi level. However, a rapid depolarization is seen with increasing temperature, 
which has been attributed to the emergence of non-quasiparticle states in the minority- 
spin band at the Fermi level.- These non-quasiparticle states offer strong evidence of the 
vital role played by electron-magnon interactions in these systems. Indeed, as was shown 
in Section VI, the electron-magnon self-energy results in a transfer of spectral weight 
in the minority-spin band near the Fermi energy. In the following, we will investigate 
correlation effects on spin stiffness of Co2MnSi within a simplified picture. 

The lattice structure of Co2MnSi is shown in Fig. |9l The Co atoms occupy the corners 
of a simple cubic lattice, and the centers of the cubes are occupied by Si and Mn atoms 
alternately. Since the sp electrons of Si do not participate in the electronic and magnetic 
properties of the system as they lie much below the Fermi level ,-22 we consider only the Co 
and Mn atoms, the hybridization of whose orbitals results in the minority-spin band gap 
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FIG. 11: Electronic DOS of Co2MnSi as obtained from LDA+DMFT studies.-^. 

at ej?.— Thus the lattice effectively consists of an admixture of the sc and bcc lattices, 
and so we consider the following simplified dispersion in our analysis: 

ek = (1 - r/)6^^^ + r/e^^ (39) 



where r/ is a tunable parameter which controls the relative strengths of the bcc and sc 
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FIG. 12: Renormalized spin stiffness corresponding to tlie band dispersion given by Eq. (|39p 
witli t = 0.15 eV and t'/t = 0.5, and the effective quantum parameter 1/A^ = 1/3. The spin 
stiffness passes through a maximum at around r] = 0.3, with a crossover to negative values 
around r/ ~ 0.7. 

characters of the dispersion, and the bcc and sc band dispersions are given by: 

e^'^^ = 8tcos{kxa)cos{kya)cos{kza) + 2t' cos(2A;^a), (40) 

= 2t{cos{2k^a) + cos{2kya) + cos{2k^a)). (41) 

The lattice parameter for Co2MnSi is 5.654 A,— nearly twice that of iron. For t = 0.1 
eV and t'/t = 0.5, Fig. [TO] shows the electronic density of states for different rj values. 
The transition from a purely bcc band for r/ = to a purely sc band for rj = 1 is marked 
by a sharp peak near the lower band edge for t] ~ 0.3, which should be particularly 
favourable for ferromagnetism. While the spectral distribution is qualitatively similar to 
the majority-spin DOS obtained from first principles calculations^ shown in Fig. [Ill the 
bandwidth (~ 1.2 eV) is underestimated by roughly 50 

We have again taken the interaction parameters U/t = 16 and J/U = 1/4 as in Section 
IV, which are similar to those considered in LSDA+DMFT studies.^ With the effective 
quantum parameter 'h' = 1/N = |, the renormalized spin stiffness as a function of band 
filling n is shown in Fig. [12] for different 1] values. Here the hopping term is taken as 
t = 0.15eV so that the majority-spin bandwidth of about 2eV matches with the first- 
principles calculations shown above.- The spin stiffness maximum at r] ~ 0.3 corresponds 
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to the sharply peaked electronic DOS in Fig. [101 and the peak value of about 600 meVA^ 
compares favourably with the stiffness value 575 meVA^ obtained^^i^i for the L2i ordered 
phase (structure shown in Fig. [9]). For the B2 phase of Co2MnSi, where the Mn and Si 
atoms are randomly disordered, a much lower spin stiffness of 324 meVA^ was obtained.— 



VI. REALISTIC MULTI-BAND SYSTEMS 

The quantum corrections analysis discussed so far was restricted to an A/'-fold degen- 
erate Hubbard model. In the following, we present an extension of our Goldstone-mode 
preserving approach to realistic multi-band models with non-degenerate orbitals and inter- 
orbital hopping, as obtained from first principles calculations. An approximate scheme 
for the evaluation of the correlation-induced quantum corrections for a general multi-band 
model is discussed below. 

We consider the following multi-band Hamiltonian for realistic systems which includes 
non-degenerate orbitals and inter-orbital hopping: 

-f^ = - ^ ti^al^aj^ - ^ U^^Si^ ■ ^^'^ = X] ^r'^icM^ki/ - ^ U^^Si^ ■ Si^ (42) 

(ij)l-iv ^jM,*^ k,/^i/ itfJ-jf 

where the generalized interaction matrix Uf^,y includes the intra-orbital Coulomb interac- 
tion Ufj, for /i = z/ and the inter-orbital Hund's exchange J^i, between orbitals fi and u for 

As in the A/'-fold orbitally degenerate case,^® the transverse spin-fluctuation propagator 
components are expressed in terms of the irreducible particle-hole propagator as: 

[x'^{^,^)]f^<^ = [<l>{^,^)]f^u + X^[0(q,c^)]MM'[^]M'<^'[x"^(q,t^)]v'v (43) 

which can equivalently be written in a matrix form in the orbital basis: 

fx-+(q 00)] = ^'^^^'^^J (44) 
^"^'^^J l-[[/][0(q,u;)]- 

Correlation-induced quantum corrections to the irreducible particle-hole propagator 
can be obtained by systematically incorporating self energy and vertex corrections as in 
the orbitally degenerate case. An approximate diagrammatic scheme for evaluating the 
first-order quantum correction 0*^^)(q, cj) in the general multi-band case is shown in Fig. 
which follows from the spin-charge coupling structure of the quantum correction and 
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FIG. 13: The spin-charge couphng structure of the first-order quantum correction (/)(^) (q, in 
the general multi-band case with non-degenerate orbitals and inter-orbital hopping. 

the magnon self energy for the saturated ferromagnet in the A/'-orbital Hubbard model.— 
The corresponding expression for the first-order quantum correction is obtained as: 

[0«(q,a;)U = 5^ /■ ^ [x;;(Q, 1^) f/^J(k; q - Q, a; - ^])] T,, (45) 

k,Q 

where the minority-spin effective interaction U^, = U^u^^y^iUu' involves the exchange of 
charge fluctuation propagator in the majority-spin bands: 

/7 / 
^[G°(k + q - Q, u;' + u; - n)]^,, [G°(k, 



if orbital (band) mixing is neglected. Furthermore, in the orbitally degenerate case, the 
effective interaction term yields the factor [[/^ + (A^ — 1) J^], and the effective quantum 
parameter is recovered as earlier. This correlation-induced coupling between the spin and 
charge fluctuations represents scattering of a magnon into intermediate spin-excitation 
states accompanied by charge fluctuations in the majority spin band. These intermediate 
states include both the sharp magnon excitations and the Stoner excitations spread over 
the Stoner continuum. Finite imaginary part of this magnon self energy due to the 
gapless charge excitations results in finite magnon damping at zero temperature even for 
low-energy magnon modes lying within the Stoner gap. 
The spin-charge interaction vertex above is given by: 



r^(k; q, Q, fi) = (^xj(k; q, w) - ^. 



Q,fi) 



(47) 
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where (in the absence of band mixing) 

(k; q, u;) = ^ j;^ , (48) 

and 

2A;(q.L;Q.n) - E4(k;q.-)4(k';Q,n)/x2(Q.(J) , (49) 

which generally has weak momentum dependence due to the averaging over momentum k'. 
For q = 0, both terms in Eq. fH7|) reduce to 1/ (2A+u)). The spin-charge interaction vertex 
r and the magnon self energy therefore vanish identically, and the Goldstone mode is thus 
explicitly preserved. For small q, F^ ~ (q.Vck)^, indicating short-range interaction. Also, 
the spin-charge coupling results in a quantum correction only to the exchange contribution 
to the spin stiffness as required; quantum corrections to the delocalization contribution 
of the type (q.V)^ek cancel exactly.— 

The RPA-level description used above in Eq. ( H5l) is obtained by replacing the irre- 
ducible particle- hole propagator [(p]^u in Eq- dlSD by the bare particle-hole propagator: 

— 5^[G?(k',c.O]M.[G5(k' - q,u;' - c.)],. (50) 
in terms of the HF-level Green's function matrices in orbital space: 

[G'^{\.M] = [^^-H^\n-\ (51) 

which are no longer diagonal due to orbital mixing by the hopping term. The HF-level 
Hamiltonian matrix in momentum-orbital-spin space: 

(k)] = [ef] - 2A, • S, (52) 

in terms of the bare band energies and the self-consistent mean fields for orbital /i: 

A^ = t/^(S^) + 5^J^,(S,) (53) 

V 

involving both the intra-orbital Hubbard interaction and the inter-orbital Hund's coupling. 
Assuming the magnetization to lie along the z direction, the exchange splitting is obtained 
as 2A^ = U^m^ + J^ui^ui where is the magnetization in orbital /x. 
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VII. CONCLUSIONS 

The correlated motion of electrons in metallic ferromagnets was investigated in terms 
of a realistic Hubbard model with Af-iold orbital degeneracy and arbitrary intra-orbital 
Coulomb interaction U and inter-orbital Hund's exchange J. A spin-rotationally sym- 
metric and non-perturbative scheme was developed to study correlation-induced quantum 
corrections beyond the RPA, wherein self energy and vertex corrections were incorporated 
systematically so that the Goldstone mode is explicitly preserved order by order. An effec- 
tive quantum parameter 'h' — ^^^Ij^Zlyy was obtained which determines, in analogy with 
1/S for quantum spin systems and 1/N ior the A^-orbital Hubbard model, the strength of 
quantum corrections to spin stiffness and magnon energies. The rapid suppression of this 
quantum parameter with Hund's coupling J, especially for large Af, provides fundamental 
insight into the phenomenon of strong stabilization of metallic ferromagnetism by orbital 
degeneracy and Hund's coupling. 

The above approach was illustrated for the case of ferromagnetic iron. Electronic and 
magnetic properties were investigated on a bcc lattice with realistic parameters. The 
electronic spectral function renormalization was seen to be especially strong near the X 
and R points, indicating significant band flattening and mass renormalization, and also 
strong non-quasi-particle character of the low-energy minority-spin states near the Fermi 
energy. With the same set of parameters, both the calculated spin stiffness and Curie 
temperature values obtained were in quantitative agreement with measurements. An im- 
portant finding of our explicitly Goldstone mode preserving scheme was the result that 
including the contribution of long wavelength modes yielded a nearly ~ 25% reduction in 
the calculated Curie temperature, which should be of interest in view of the significant 
over-estimation of the Curie temperature in approaches where only local magnetic exci- 
tations are included. Another illustration was provided for the half metallic Heusler alloy 
Co2MnSi and the calculated spin stiffness was again in agreement with measured values. 

Finally, an outline was presented for calculating the correlation-induced quantum cor- 
rections for a generic multi-band metallic ferromagnet including realistic band-structure 
features of non-degenerate orbitals and inter-orbital hopping as obtained from LDA cal- 
culations. This extension drew upon the spin-charge coupling structure, which explicitly 
preserves the Goldstone mode, highlights the role of majority-spin charge fluctuations on 
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spin dynamics, and also explicitly yields finite zero-temperature magnon damping due to 
decay into longer wavelength modes accompanied with internal charge excitations. This 
provides an important advantage over other approaches which incorporate correlation 
effects in terms of an effective Heisenberg model only. 



VIII. APPENDIX 



The identity used in Eq. fl20p is derived here. From Eqs. fll8fl9p and with U 
[U+ + (A/" - 1)U-]/^ and J = [f/+ - U-]/^, we obtain: 



1 



t/r;+ + (Ar-i)jr;+ = -^ 



1 - u+xo 1 - u-xo 



which can be written in terms of a superposition: 

1 



1 
If 



1 - f/+Xo 1 - U-x, 



- A 

~ M 
B 

+ AT 



A/--! 



1-U+xo l-U-xo, 
1 1 



1-U+Xo l-U-Xo. 
of the two functions in Eqs. (|T2l) and ( IT3l) . Solving for A and B yields: 



B = 2(A/'- l)f/J+ (A/"- l)(A/'-2)j2, 

so that in terms of Xaa Xap Eqs. ffT^ and ffT^ . we obtain: 

UT-^ + (A/" - i)JKt = {[U' + {M- i)Ax-i 

+ [2{U - l)UJ 2)J^]x-;}/x,. 
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